DIQUARK CONDENSATION IN TWO COLOUR QCD 
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Unqucnched lattice SU(2) is studied at nonzero chemical potential in the strong 
coupling limit. The topic of diquark condensation is addressed analyzing the prob- 
ability distribution function of the diquark condensate. We present results at zero 
external source without using any potentially dangerous extrapolation procedure. 
We find strong evidences for a (high density) second order phase transition where 
a diquark condensate appears, and show quantitative agreement of lattice calcula- 
tions with low-energy effective Lagrangian calculations. 

1 Introduction 

Recently the standard scenario for the phase diagram of QCD in the chemi- 
cal potential-temperature plane has changed; besides the standard hadronic 
and quark-gluon plasma phases., the existence of a new state of matter has 
been claimed by several groupscl. This new phase is characteristic of the high 
density-low temperature regime; the asymptotic freedom of QCD and the 
known instability of large Fermi spheres in presence of (whatever weak) at- 
tractive forces results in a pairing of quarks of higher momenta (in analogy 
with the Cooper pairing in solid state systems at low temperature). A con- 
densation of quark pairs should be the distinctive signal of the new phase and 
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has been indeed predicted using simplified phenomenological models of QCD. 

Unfortunately the lattice approach, the most powerful tool to perform 
first principles non perturbative studies, is affected in the case of finite density 
QCD by the well known sign problem that has prevented until now any step 
towards the understanding of this new physics. This is not the case for the 
SU(2) theory where the fermions are in the pseudo-real representation and 
finite density numerical simulations are feasible. 

In the following we present a detailed study of diquark condensation for 
the unquenched two colours model. This work follows a previous one where 
we considered chiral and diquark susceptibilities and found strong evidence 
for a phase transition separating the ordinary low density phase from a high 
density one where chiral symmetry is restored and the baryon number sym- 
metry is broke no. In the present work we mainly focus on an approach, based 
on the analysis of the probability distribution function (p.d.f.) of the order 
parameter, to extract the value of the order parameter directly at zero exter- 
nal source without using any potentially dangerous extrapolation procedure, 
implicit in the standard approach. 

The study of diquark condensation in two colours QCD is an interesting 
topic by itself and, notwithstanding the relevant differences between SU(2) 
and SU(3) at non zero baryon density, we can hope to use some of the results 
to get insights about the three colours case. Even if the phase diagram of 
SU(3) and SU(2) are expected to share similarities, we have to be aware of 
the differences between the two models. 

• The quark-quark condensate {qq} is coloured for the N c = 3 theory and 

colourless for the N c — 2 one. This implies that in the SU(3) case the 
diquark condensation has to be interpreted as an Higgs-like mechanism 
leading to the breaking of the local colour symmetry, while in the SU(2) 
case we have a standard spontaneous symmetry breaking (SSB) of the 
Uv(X) symmetry associated with conservation of baryon number. 

• The zero temperature critical chemical potential, related to the mass of 

the lightest baryonic state, is different in the two cases: it is 1/3 of the 
nucleon mass for SU(3) and 1/2 of the pion mass for SU(2) (i.e. zero at 
vanishing quark mass). For N c = 2, in order to have a phase diagram 
similar to the SU(3) one, it is appropriate to consider a non zero bare 
quark mass. 

• The lightest baryonic state is a fermion for the tree colours model and a 

boson for the two colours one. 

Even if our results can not be directly extended to the SU(3) case we can 
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use lattice calculations to do a highly non trivial, quantitative check of the 
continuum predictions for the SU(2) theory. In particular we have considered 
results coming from low energy effective Lagrangian calculations! 

In the next section we present some details of the algorithm. The last 
section is devoted to the presentation and discussion of numerical results. 

2 Overview of the numerical technique 

2.1 Simulation scheme 

The standard way to study SSB is to introduce first an explicit symmetry 
breaking term in the action. If we do that for the diquark in the SU(2) model 
we have to add a terma j(ipT2ip + 4>T2'4>) and, after integrating the Grass- 
mann field, the fermionic contribution for Nf = 4 quark flavours becomes 
proportional to the Pfaffian of a 4V x 4V matrixH: 



and A is the usual lattice Dirac operator (it contains the mass and [i de- 
pendence). Using the relation T2AT2 = A* we can easily prove that B is 
antihcrmitian and det(B + j) > for any j. 

It can also be shown that the eigenvalues of B are doubly degenerate and 
B 2 is in a block diagonal form with two hermitian blocks on the diagonal 
having the same eigenvalues. It follows that to compute det(£? + j) for any 
j (i.e. to obtain all the eigenvalues of B) it is sufficient to diagonalize only 
one block of B 2 (reducing the problem to the diagonalization of a 2V x 2V 
hermitian matrix) and then take the two pure imaginary square roots of the 
(real and negative) eigenvalues. 

To avoid the sign ambiguity in (Q) it is customary to consider a theory 
with Nf — 8 quark flavours where the fermionic partition function becomes 
Zf e rm{j) = det(B + j). On the other side if we are able to work directly at 
zero diquark source things become much simpler. In the j = limit the sign 
ambiguity disappears and the Pfaffian is positive definite since Pf(B) = det A 
and the last quantity is real and positive for any value of \i. Then we can 
easily consider any value of Nf writing Zf erm (j — 0) = (det(B)) N f / 8 . In the 
next section we show how it is possible to take advantage of this feature and 
extract the value of the order parameter for the diquark condensation for all 
Nf without any extrapolation to j = 0. 
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Z ferm (j)=Pf(B + j) 



±^det(B + j) 



(1) 



where B is 





We have studied the phase structure of the theory in the limit of infinite 
gauge coupling ((3 — 0). To simulate the (3 = limit of the theory we have 
measured fermionic observables on gauge configurations generated randomly, 
i.e. with only the Haar measure of the gauge group as a weight. This choice 
implies a Gaussian distribution of the plaquette energy around zero which, 
according to the results of Morrison and Handsel, has a net overlap with the 
importance sample of gauge configurations at the values of /x and m used in our 
calculations. The validity of this procedure has also been tested comparing 
different physical observables (number density and chiral condensate) with 
Hybrid Montecarlo resultsO. 

We have considered the theory in a 4 4 and 6 4 lattice diagonalizing 300 
gauge configurations in the first lattice volume and 100 in the second one. 
As we pointed out in the introduction, to stay closer to the SU(3) case the 
simulations have been performed at non zero quark mass. We choosed m = 
0.025,0.05,0.20 and values of the chemical potential ranging from /i = to 
fx = 1.0. 

All numerical simulations have been performed on a cluster of Pentium 
II and PentiumPro at the INFN Gran Sasso National Laboratory. 

2.2 Analysis of the probability distribution function 

The use of the p.d.f. to analyze the spontaneous symmetry breaking in spin 
system or Quantum Field Theories with bosonic degrees of freedom is a stan- 
dard procedure. Less standard is its application to QFT with Grassmann 
fields where, for obvious reasons, the fermionic degrees of freedom have to 
be integrated analytically. Nevertheless this method has been developed to 
extract the chiral condensate in the chiral limit from simulations of QFT with 



The same ideas can be used to study the vacuum structure of two colour 
QCD at non zero density and specifically to extract the diquark condensate 
at j = 0. We refer to the original paper for a full description of the p.d.f. 
technique and present a brief introduction focusing on the peculiarities of the 
diquark condensate case. 

Let a be an index which characterizes all possible (degenerate) vacuum 
states and w a the probability to get the vacuum state a when choosing ran- 
domly an equilibrium state. If c a is the order parameter in the a state we 
can write 



where V is the lattice volume and the sum is over all lattice points. P(c), the 
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p.d.f. of the diquark order parameter c, will be given by 

P ( C ) = ^WqS(C- Cg) 

a 

= lim — f[dU][d^][di;}e- SG ^ + ^S{^y^ ^t 2 ^(x) + ^pT 2 i;(x) - c) 

X 

The main point is that while P(c) is not directly accessible with a numer- 
ical simulation its Fourier transform 

P(q) = J dce iqc P(c) (4) 

can be easily computed. Inserting in (Q) the definition of P(c) and using an 
integral representation for the 5-function, after some algebra we arrive at the 
following formula 

where Pv{l) is the Fourier transformed p.d.f. of the diquark order parameter 
at zero external source and finite volume. 

To take advantage of (g) we should be able to compute correctly the 
Pfafnans involved. This indeed turns out to be easy. First note that, as 
previously stated, PfB(j = 0) is a real and positive quantity so that the only 
ambiguity is related to the sign of PfB(j = y) = ±y/dctB(j = iq/V). If we 
remember that the 4V eigenvalues of the antihermitian matrix B are doubly 
degenerate and can be written as ±zA„ (n — 1, • • • , V) with real and positive 
A„ we arrive at the following expression 

n— 1 

The sign ambiguity can now be solved noticing that has to be positive at 
j = and, increasing q, changes sign each time j is equal to one of the eigen- 
values. The only possible exception would be an (accidental) extra degeneracy 
in the eigenvalues of B resulting in the Pfafnan not crossing zero but tangent 
to the horizontal axis. This situation never occurred in our simulations. 

Once we have the Pv(q) and hence, by simple Fourier transform, the 
p.d.f. of the order parameter we need to extract the correct value for the 
order parameter. To do that we first have to recall that diquark condensate is 
the order parameter of the Uy(X) symmetry associated with baryon number 
conservation. Indeed (V'^V' + V ,T 2'0) an d i{ipT2ip — ipT2ip) are the components 
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Figure 1. P(c) for the 4 4 lattice at fi = 0.2,0.4,0.475,0.5 (from top to bottom). 

of a vector (in a plane) which rotates by an angle 29 when we do a global 
phase transformation of parameter 9 on the fcrmionic fields. Therefore, if cq 
is the vacuum expectation value of the diquark condensate corresponding to 
the a- vacuum selected when switching-on a diquark source term, P(c) can be 
computed as 



which gives P(c) — 1/(7t(cq — c 2 ) 1 / 2 ) for — c < c < c and P(c) = 
otherwisea. In the symmetric phase cq = and P(c) reduces correctly to 
a (5-function in the origin. 

The above results are valid in the thermodynamic limit while, at finite 
volume, the non analyticities of the p.d.f. are absent. Without entering in the 
details of a finite size scaling analysis we expect, for the finite volume p.d.f. 
Py(c), a function peaked in the origin in the symmetric phase and peaked at 
some non zero value in the broken phase. This is indeed the behaviour we 
can observe in fig. [j] where the the p.d.f. of the smallest volume is reported 
at different values of fi. It is clear that, increasing the chemical potential, 
the vacuum starts to be degenerate signalling a spontaneous breaking of the 
baryon number conservation. 

We can also compare the -Py(c) for two lattice volumes in the symmetric 
and broken phase. This is done in fig. 2 where we see clearly as, increasing 
the volume, the peak of the p.d.f. becomes sharper. To determine the value 
of the diquark condensate we used the position of the peak: a definition that 
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clearly converges to the correct value in the thermodynamic limit. 

From fig. (||) we see also that data of the larger volume are more noisy 
(indeed we also get negative values for the p.d.f. in the broken phase). To 
have an estimate of the errors we calculated several distribution functions for 
the 6 4 lattice using independent subsets of our data. We saw that the position 
of the peak was very stable, within 1 percent in the broken phase, and we used 
this quantity for the errorbars. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.0 0.2 0.4 0.6 0.8 1.0 1.2 



Figure 2. P.d.f. for 4 4 (dashed line) and 6 4 (continuous line) in the symmetric (left) and 
broken (right) phase. 



3 Results 

Here we present the results for the diquark condensate at j = as a function 
of the chemical potential. Fig. 3 contains our data points for the 6 volume 
for several values of Nf at two different masses: m = 0.025 and m = 0.2 . 

First we can derive a peculiar phase diagram: two symmetric phases 
separated by a broken one and two possibly continuous transition points. This 
result is not new and has already .been predicted by Mean Field calculations^] 
as well as numerical calculationso. The high density symmetric phase has 
no physical relevance, it is simply the consequence of the saturation of all 
lattice sites with quarks. This saturation phenomena has nothing to do with 
continuum physics and is a pure lattice artifact. The physically interesting 
phase transition, i.e. the transition that has a continuum counterpart, is the 
first one. 

To test our numerical results we have considered the analytical predictions 
of low energy effective Lagrangian. The j ' = case is described by a simple 
formula where the only two parameters are the chiral condensate and one half 
the pion mass calculated at j = 0, jjb = Od. We did independent simulations 
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Figure 3. Diquark condensate as a function of fi for 6 4 lattice (points) and the continuum 
predictions (line) for m = 0.025 (left) and m = 0.2 (right). 



Table 1. Parameters for the low energy effective Lagrangian predictions. 



m 




m,/2 


0.025 


1.31(2). 


0.1696(11) 


0.2 


1.23(2) 


0.4841(7) 



to calculate this quantities at the same values of m used in the previous 
analysis, obtaining the values in Table 1. Another prediction of the continuum 
calculation is that the value of the diquark condensate is independent of the 
flavour number. 

The continuous line in fig. 3 is the analytical prediction for this model. 
We see a remarkable agreement with our strong coupling results up to values 
of the chemical potential where the saturation effects start to be relevant. We 
see also that the Nf dependence in the simulation points is small. 

We checked explicitly the low energy effective Lagrangian prediction at 
finite j too. In that case we used standard techniques to compute chiral and 
diquark condensate as well as the number densitycl in a Nf — 8 simulation 
(fig. 4). Also in that case (where no phase transition is present and the 
predictions are still independent of Nt ) the agreement is very good up to 
values of \l where the number density becomes a consistent fraction of its 
maximum lattice value (i.e. (J>/m n ~ 0.6). 

It is surprising that a (3 = calculation has a so good agreement with a 
continuum prediction. Since the analytical predictions are, at the values of 
\i presented, well inside the validity region of the low energy approximation 
(fi << the mass of the first non goldstone excitation) we can conclude that 
this gives an indication of a very small dependence of lattice results on (3. 
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Figure 4. Low energy Lagrangian predictions (lines) and our Nf = 8 data (symbols) for 
chiral condensate, diquark condensate and number density at m = 0.2 and j = 0.1m vs. 
Ii/m-jr- 6 4 data are squares and circles, 4 4 data are diamonds. 
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